Abstract. This is the second paper in a series of three 
in which we present an exhaustive inventory of the sohd 
state emission bands observed in a sample of 17 oxygen- 
rich dust shells surrounding evolved stars. The data were 
taken with the Short and Long Wavelength Spectrographs 
on board of the Infrared Space Observatory (ISO) and 
cover the 2 to 200 /im wavelength range. Apart from the 
broad 10 and 18 /zm bands that can be attributed to amor- 
phous silicates, at least 49 narrow bands are found whose 
position and width indicate they can be attributed to crys- 

. talline silicates. Most of these emission bands are concen- 

I trated in well defined spectral regions (called complexes). 

i We define 7 of these complexes; the 10, 18, 23, 28, 33, 
40 and 60 micron complex. We derive average properties 
^ of the individual bands. Almost all of these bands were 
K-^ ' not known before ISO. Comparison with laboratory data 

. suggests that both olivines (Mg2a;Fe(2_2a;)Si04) and py- 
,—1 ' roxenes (Mgj;Fe(i_a;)Si03) are present, with x close to 1, 

■ i.e. the minerals are very Mg-rich and Fe-poor. This com- 

■ position is similar to that seen in disks surrounding young 
^ stars and in the solar system comet Hale-Bopp. A signif- 

' leant fraction of the emission bands cannot be identified 
, with either olivines or pyroxenes. Possible other materi- 
T— I ' als that may be the carriers of these unidentified bands 
, are briefly discussed. There is a natural division into ob- 

■ jects that show a disk-like geometry (strong crystalline 
-«.^ ■ silicate bands), and objects whose dust shell is charac- 

' teristic of an outflow (weak crystalline silicate bands). In 
particular, stars with the 33.5 /zm olivine band stronger 



Oh- 



O ' than about 25 percent over continuum are invariably disk 
^ , sources. Likewise, the 60 fim region is dominated by crys- 
^ talline silicates in the disk sources, while it is dominated 
• • : by crystalline H2O ice in the outflow sources. We show 
^ ' that the disk and outflow sources have significant differ- 
, ences in the shape of the emission bands. This difference 

must be related to the composition or grain shapes of the 
5^ , dust particles. The incredible richness of the crystalline 

silicate spectra observed by ISO allows detailed studies of 

the mineralogy of these dust shells, and is the origin and 

history of the dust. 
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1. Introduction 

At the end of their life both low and high mass stars loose 
a large fraction of their mass in the form of a dense stel- 
lar wind. When the temperature in the outer regions of 
the atmosphere becomes low enough solid material {dust) 
condenses out of the gas. Mass loss can eventually dom- 
inate stellar evolution. Since the dust may also play an 
important role in the mass loss process, it is interesting to 
investigate the physical and chemical processes responsi- 
ble for dust formation. One way to do this is to study the 
endproducts of this dust formation process. The compo- 
sition of the dust which has condensed provides valuable 
information on the conditions when the dust was formed 
and thereafter. 

The dust around stars can be observed at different 
wavelengths. In the visible and near-infrared (NIR) one 
can look at the wavelength dependence of the absorption 
caused by the dust. Another way to investigate the dust 
at these wavelengths is by means of scattered light of the 
central source by small dust particles in the circumstel- 
lar environment or with polarimetric observations. In this 
paper we will study the emission of the dust in the mid- 
infrared (MIR) and far-infrared (FIR). 

The dusty environments around evolved stars can be 
divided in carbon-rich (C-rich) and oxygen-rich (0-rich) 
environments, depending on the C/0 ratio of the mass- 
losing stars. This division is the result of the stability of 
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the CO molecule, which is formed before the dust con- 
denses. If there is more carbon than oxygen (C/0 > 1) all 
the oxygen will be trapped in CO, and the dust species will 
be carbonaceous, e.g. SiC, Polycyclic Aromatic Hydrocar- 
bons (PAHs) or amorphous carbon. If the C/0 ratio is 
smaller than one, i.e. there is more oxygen than carbon, 
all the carbon is trapped in CO and 0-rich dust will be 
formed, e.g. simple oxides and silicates. 

In 1995 the Infrared Space Observatory (ISO; Kessler 

et al. 1996) was launched which opened the possibility to 
study mass-losing stars at infrared wavelengths with un- 
precedented wavelength coverage and spectral resolution. 
Before the launch of ISO it was generally assumed that in 
the dusty winds of 0-rich evolved stars only amorphous 
silicates were formed. One of the remarkable discoveries 
of ISO was the detection of crystalline silicates outside 
our own Solar System. Their spectral signature was not 
only found in the spectra of young stars (Waelkens et al, 
1996), but also in the outflows of evolved stars (Waters 
et al, 1996). In the latter case, we even found one exam- 
ple, IRAS09425-6040, where 75% of the circumstellar dust 
consists of crystalline silicates (Molster et al. 2001). 

In contrast to amorphous silicates, crystalline silicates 
provide a unique opportunity to determine for the first 
time the chemical composition of the dust particles. The 
relatively sharp features of the crystalline silicates are very 
sensitive to compositional changes, this in contrast to the 
broad and smooth amorphous silicate features. Crystalline 
silicates will also help us to better understand the physi- 
cal and chemical conditions under which dust is formed. 
In particularly, they need high temperatures and a slow 
cooling down to form. 

Crystalline silicates are found in the outflows of 
evolved stars which replenish the ISM, and also around 
young stars which form from the ISM (e.g. Waters et al. 
1996; Waelkens et al. 1996). However, up to now no crys- 
talline silicates have been found in the ISM. The abun- 
dance in the ISM might be too low. However, this would 
make it difficult to explain the high abundance of crys- 
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talline silicates in a young star as HD100546 (Waelkens 
et al. 1996; Malfait et al. 1998). On the other hand, if 
the crystalline silicates are destroyed in the ISM then the 
question arises how they are formed around young stars. 
Possible solutions to this last problem are given by Molster 
et al. (1999a; hereafter MYW). 

The crystallization process is not well understood and 
in order to better constrain the origin of these grains, it is 
important to accurately describe the observed properties 
in different environments. After the first reports of the 
crystalline silicates, only a few objects have been anal- 
ysed in some detail, e.g. HD100546 (Malfait et al. 1998), 
CPD-56°8032 (Cohen et al. 1999), AFGL4106 (Molster 
et al. 1999b), IRAS09425-6040 (Molster et al. 2001). This 
will be the first attempt to obtain an average spectrum of 
the different crystalline silicate features which are found 
around evolved stars. 

In the first paper of this series (Molster et al. sub- 
mitted; hereafter paper I) we have described the infrared 
spectra of 17 sources. Here, we will derive mean spectra 
based on these sources and identify the different features. 
In the third paper of this series (Molster et al. submit- 
ted; hereafter paper III) we will apply a simple dust emis- 
sion model to the spectra of our sample and derive several 
(spectral) trends. 

In Section we will briefly summarise how we obtained 
these mean spectra. For a more extensive description we 
refer to paper I. The results will be shown in Section ^, 
where we will also order and define the different solid state 
bands and complexes found. In Section |4| we will identify 
most of the features found in Section ^ by comparison to 
laboratory spectra of cosmic dust analogues. 

2. The derivation of the mean spectra. 

2.1. The sample 

We have selected 17 evolved stars which all showed the 
presence of crystalline silicates in their infrared spectra. 
These 17 sources covered a wide range in evolutionary 
status. The list included 2 AGB stars (OH26.5-h0.6 and 
IRAS 09425-6040), 5 post- AGB stars (HD44179, 89 Her, 
AC Her, Roberts 22 and HD161796), 3 planetary neb- 
ula (NGC6537, NGC6302 and VY 2-2), three massive 
stars (IRC-l-10420, AFGL4106 and NML Cyg) and 4 ob- 
jects with poorly known evolutionary status (MWC922, 
MWC300, HD45677 and HD179821). An extensive dis- 
cussion of the individual stars and their spectra can be 
found in Paper I (Molster et al. 2002a). 

For all stars in this sample, spectra were taken with the 
Short Wavelength Spectrometer (SWS; de Graauw et al, 
1996) on board ISO and for most of them we also got spec- 
tra taken by the Long Wavelength Spectrometer (LWS; 
Clegg et al. 1996). 

As we will show later in this paper, there is a clear 
difference in the spectra of stars which (are expected to) 



have a disk and those who have not. This difference is seen 
not only in the strength of the features (MYW), but as 
we will demonstrate also in the shape and position of the 
features (see below and Molster et al. 2002b, hereafter Pa- 
per HI). Therefore, we have decided to divide our sample 
accordingly into 2 groups. 10 stars in our sample belong 
to the 'disk' sources and 7 stars belong to the 'non-disk' 
or 'outfiow' sources. 

Although the spectral shape of the different complexes 
is different for the disk and the outflow sources, the indi- 
vidual bands are present in both type of classes (although 
in different strength ratios!). The average FWHM (full 
width at half maximum) and position of individual bands 
is therefore determined based on the complete sample and 
no discrimination has been made between disk and outflow 
sources. 

2.2. Mean complex spectra 

For all stars we defined a continuum, which we subtracted 
from the data to obtain a better view of the dust features; 
see Paper I for more details on how we derived the con- 
tinuum. The resulting continuum subtracted spectra were 
used to determine the mean complex spectra. 

In the (continuum subtracted) spectra of the stars, 
we can identify 7 different complexes in which multiple 
components can be identified; some of which are severely 
blended. These complexes are found near 10, 18, 23, 28, 
33, 40 and 60 /im, and named after their position, resp. 
the 10, 18, 23, 28, 33, 40 and 60 micron complexes. 

In order to get a reference spectrum for the different 
complexes we have derived 'mean' spectra for each com- 
plex. The mean complex spectra were derived for both the 
outflow and the disk-sources separately. In order to obtain 
these mean complex spectra we have extracted each com- 
plex from all the continuum subtracted spectra. The com- 
plexes of the stars belonging to one group (disk or outflow 
sources) were added, using a weighing factor proportional 
to the S/N of that part of the spectrum, to create a 'mean' 
spectrum (see Paper I for more details). 

3. Results for the complexes 

In this section we will systematically discuss the general 
trends in the spectra of our program stars, by sub-dividing 
the sample into disk-sources and outflow-sources, and by 
dividing the spectra into the different complexes: the 10, 
18, 23, 28, 33, 40 and 60 micron complex. Because of the 
low dust temperature of most sources and therefore low 
flux levels below 7 /xm and the absence of crystalline sili- 
cate features in this wavelength range, we decided to limit 
our investigation to wavelengths longwards of 7 //m. 

First, we will give an overview of the results that apply 
to the whole sample. Then, we will discuss the different 
complexes individually, after which we will identify many 
of the features. In the discussion of the different complexes 
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Table 1. The characteristics of the different features. The mean A (FWIiM) is the average wavelength (Full Width 
Half Max) of the feature, with the errors taken into account. A (FWHM) min and max are the minimum and maximum 
value found in the sample. The N indicates not a member of one of the 7 complexes 
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Table 2. The presence of the features which are seen in at least three different SWS spectra of our sample, + means 
emission, - is absorption, o means no good data available, 2 implies that a feature split in 2 separate features in this 
source and underlined x's mean all present, but measured as one feature. For MWC922 the blend of the 9.1 and 9.5 and 
for OH26.5+0.6 the blend of the 23.0 and 23.7 micron features are in absorption, while all other blends are in emission. 
The order of the stars is in decreasing strength of the blend of the 33.6 micron feature, am = amorphous silicate, d= 
diopside; (o/c)e = (ortho/clino)enstatite; f = forsterite; i = instrumental artifact; s= silica; si = a crystalline silicate; 
w = crystalline water-ice. 





CO 
00 


CO 


oj 




CS 


0) 


,4 f,d? 


0) 
CN 


O 


CN 


o 


0) 




f,e 


03 


LO 


,7 s?,d? 


LO 




o 




03 

00 


LO 


o 


1 f,s? 


00 


Bandname 


(35 


CO 

d> 


o 


o 


T— 1 


LO 
1 — 1 


lO 
1 — 1 


1 — 1 


CO 

1 — 1 


1 — 1 


00 
1 — 1 


00 
1— I 


00 


03 


o 

CN 


CN 


CN 
CN 


CO 

CN 


CO 

CN 


CO 

CN 


CN 


lO 
CCI 


CO 

CM 


CO 


IRAS09425 












? 




+ 


+ 


+ 


+ 






+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


+ 


+ 


+ 


NGC6537 






















+ 






+ 


+ 


+ 


+ 


+ 




+ 


+ 




+ 


+ 


+ 


+ 


NGC6302 




























+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


+ 






MWC922 




X 


X 










+ 






2 


+ 




+ 


+ 


+ 


+ 


+ 




+ 


+ 


+ 


+ 


+ 


+ 


+ 


AC Her 


-X- 
1 


-X- 
1 


-X- 
1 


-X- 
1 


-X- 
1 


-X- 
1 


+ 


-X- 
1 


-X- 
1 


-X- 
1 


+ 


-X- 
1 


? 


+ 


_l_ 


_l_ 


+ 


+ 


_l_ 


X 




_l_ 


_l_ 


_l_ 






HD45677 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


+ 


+ 


+ 


89 Her 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


+ 


+ 


+ 








+ 


+ 




X 


X 






MWC300 






















+ 


+ 




+ 


+ 


+ 


+ 




+ 


+ 


+ 




+ 


+ 


+ 




VY 2-2 










+ 
























+ 


+ 




X 


X 


X 


X 


X 






HD44179 


















+ 


+ 


+ 




? 


+ 


+ 


+ 


+ 


+ 




+ 


+ 




+ 


+ 


+ 




HD161796 










+ 






+ 


+ 


+ 






+ 


+ 


+ 


+ 


+ 


+ 




+ 


+ 




+ 


+ 


+ 




OH26.5+0.6 




























- 






- 






X 


X 












Roberts 22 


















X 


X 


? 


? 


+ 


? 


? 


? 


+ 


+ 




X 


X 


X 


X 


X 




o 


HD179821 










+ 








X 


X 


+ 




+ 


X 


X 


+ 


+ 




+ 


+ 


+ 




X 


X 


+ 




AFGL4106 










+ 








X 


X 


+ 




+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


+ 


+ 




NML Cyg 
















+ 


+ 


+ 


+ 




? 


+ 


+ 




+ 


+ 




+ 


+ 




+ 









IRC+10420 




+ 


+ 




+ 


? 


+ 


+ 


+ 


+ 


+ 




+ 


X 


X 


X 


+ 


+ 




+ 


+ 




+ 




+ 


























e,d? 


0) 


o 




c^- 
T3 






w,ce 


CP 

Q 


ce,d? 
















CD 


CM 


00 


CO 


CO 


CN 


CN 


00 


Oj 


CO 


T 1 


o 


o 


lO 


00 


LO 


00 


o 


00 






CO 






o 




Bandname 




00 

CM 


00 

CM 


03 
CN 


o 
CO 


1 — 1 

CO 


CO 


CN 

CO 


CN 

CO 


CO 
CO 


CO 


CO 


m 
CO 


CO 


03 

CO 


o 




CO 


CO 






00 




CO 




I— 1 
Oi 


IRAS09425 


+ 


+ 




+ 


+ 




+ 


+ 




+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 














NGC6537 


X 


X 


+ 


+ 


+ 




+ 


+ 




+ 


+ 


X 


X 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


NGC6302 


+ 


+ 




+ 


+ 


+ 


+ 


+ 




+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


MWC922 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


X 


X 


+ 


+ 




+ 


+ 


+ 


+ 


+ 


X 


X 


+ 


+ 


+ 




AG Her 


+ 


+ 




+ 


+ 


7 


+ 


+ 




X 


X 




X 


X 


X 


X 






+ 
















HD45677 


+ 


+ 




+ 


+ 




+ 


+ 




+ 


X 


X 


X 


X 


+ 


+ 


+ 


+ 


+ 


+ 














89 Her 


+ 


+ 




? 








+ 




+ 


X 


X 






























MWC300 


X 


X 


+ 


+ 


+ 


+ 




+ 




+ 


+ 


X 


X 


X 




+ 




X 


X 


X 














VY 2-2 








+ 


+ 






+ 




+ 


+ 


+ 


+ 


+ 




+ 


X 


X 


X 








+ 


+ 




+ 


HD44179 


+ 


+ 




+ 


+ 




+ 


+ 




+ 


+ 


+ 


+ 


+ 




+ 






+ 




X 


X 


+ 


+ 


+ 




HD161796 


+ 


+ 


? 


+ 


+ 




+ 


+ 




+ 


+ 


+ 


X 


X 


+ 


+ 


+ 


+ 


? 




X 


X 


+ 


2 


+ 




OH26. 5+0.6 








+ 








X 


X 


X 


X 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


X 


X 


+ 


+ 


+ 




Roberts 22 


+ 


+ 


+ 











+ 




+ 


+ 


+ 


X 


X 


X 


X 


+ 


+ 


+ 










+ 


7 


+ 


HD179821 


+ 


+ 




+ 






+ 


+ 




+ 


+ 


+ 


+ 




+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 




AFGL4106 


+ 


+ 












+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


X 


X 




+ 






NML Cyg 














+ 




+ 


+ 


+ 


+ 


+ 


+ 


+ 




X 


X 


+ 


+ 


+ 


+ 


X 


X 




+ 


+ 




IRC+10420 


+ 


+ 













+ 


+ 


+ 




+ 


+ 




? 


+ 


+ 


+ 




+ 


+ 


+ 


+ 


+ 







6 



F.J. Molster et al.: The crystalline silicate complexes 



and features, we will refer to feature names as xx micron 
features and to wavelength positions as yy /im, e.g. the 
30.6 micron feature of IRAS09425-6040 is located at 30.48 
fiia (Paper I). 

In Table 0, we have grouped all features seen in at 
least 2 stars and indicated their wavelength and FWHM 
limits. In principle, we assigned one feature to one spec- 
tral structure. If in at least two stars a spectral structure 
was split into two different components, we assigned two 
features to it. If in other stars these two features could not 
be separated we assumed them to be heavily blended and 
measured them together. The one exception is the 16.9 
micron feature. In MWC 922 it was found to be split into 
two different components, while in all other stars it was 
blended. We allowed the features to show a small shift in 
wavelength position. The ranges in shifts are indicated in 
Table |l| by A,nin and Amax- This method leads to a mini- 
mum number of features to analyse all spectra. Measured 
blends and uncertain features were not taken into account 
for the values in Table |l|. 

In the same table we also report to which mineral the 
feature could be traced; This will be discussed extensively 
in Section ^. We decided not to discuss the C-rich features 
which are found in many of the disk sources. 

The mean spectra of the different complexes are shown 
in Fig. 0to Fig. 0. In the next subsections we discuss these 
different complexes for all stars. We point out that the in- 
dividual disk source and outflow source spectra show sig- 
nificant source to source variations in shape and strength 
of the spectrum. Nevertheless, clear trends can be seen 
that are most obvious in the mean spectra. 

Table || gives a schematic quick-look overview of what 
features are found in what star. We indicate by using dif- 
ferent symbols the blends and the features which are un- 
certain. The order of the stars is in decreasing strength of 
the blend of the 33.6 micron feature. 

3.1. 10 micron complex (7-13 fim) 

In Fig. the mean 10 micron complexes are shown. The 
mean disk spectrum is only based on the 3 stars which 
show an 0-rich dust spectrum in this area, AC Her, 
IID45677 and 89 Her. All other disk sources have contribu- 
tions of carbon-rich material. The mean outflow spectrum 
is based on the outflow sources which show the silicates in 
emission, i.e. we have not used the absorption spectra of 
NML Cyg and OH26.5-I-0.6, because they have too strong 
absorption in this wavelength range. The mean outflow 
spectrum is dominated by a broad feature, that peaks at 
about 10.6 /im. 

In the mean disk spectrum, several narrow features are 
found, at 8.3, 9.14, 9.45, 9.8, 10.7, 11.05 and 11.4 /im. We 
note that it is likely that there is an underlying broad band 
in the 10 micron complex spectra of the disk sources. Its 
shape and peak position seem somewhat different from 
that seen in the outflow sources. However, we want to 




8 10 12 8 10 12 

"wavelength (fJ-m) 



Fig. 1. The normalised mean 10 micron complex spec- 
trum of the disk and outflow sources. The normalised 
mean spectrum for the disk sources is derived from AC 
Her, HD45677 and 89 Her, the normalised mean spectrum 
for the outflow sources is derived from all spectra except 
OH26. 5-1-0.6 and NML Cyg. 




J n- J 
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wavelength (/im) 

Fig. 2. The normalised mean 18 micron complex spectrum 
of the disk and outflow sources. All stars, except NML Cyg 
and OH26.5-f0.6 for the outflow sources, are used. 

stress that significant source to source variations are ap- 
parent in the individual spectra (see Paper I). 

3.2. 18 micron complex (15-20 ^m) 

In Fig. 1^, the mean 18 micron complexes are shown for the 
disk and outflow sources. OH26.5-I-0.6 has both the amor- 
phous and crystalline silicates in absorption in the 18 mi- 
cron complex. Therefore, we have excluded this star from 
the mean spectrum. Also NML Cyg shows the amorphous 
silicate in absorption in the 18 micron complex, although 
the crystalline silicates seem already in emission we have 
neglected this star for the mean outflow spectrum too. Ex- 
cluding the two high flux, high S/N sources IRC+10420 
and AFGL4106, did not modify the average significantly. 
We note that some residual fringing is still visible in the 
average profile. 

The 18 micron complex can be divided into narrow 
features at 15.2, 15.9 (with a shoulder at 16.2), 16.5, 16.9, 
17.5, 18.0, 18.9 and 19.5 /im. In both classes, the same nar- 
row features are found. In all sources, there are at least 
indications for the presence of amorphous silicates. How- 
ever, the outflow sources show a smoother complex due 
to the larger abundance of amorphous silicates (see also 
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Fig. 3. The normalised mean 23 micron complex spec- 
trum of the disk and outflow sources. All sources are used 
except for OH26.5-I-0.6, which still has all the features in 
absorption. 

MYW). The presence of the crystalline silicates (the nar- 
row peaks) are most prominent in the disk-sources. Note 
the weakness or absence of the 19.5 micron feature in the 
outflow sources. 

The 18.1 micron feature might be a blend of 3 features. 
There is an indication for this in the spectrum of HD44179. 
The triple structure found in this spectrum is seen in all 
datasets available for this object. Also the mean spectrum 
has a hint of these three features, indicating that they 
might be weakly present in the individual spectra, close 
to the detection limit. 

3.3. 23 micron complex (22. - 25.5 fim) 

The mean 23 micron disk and outflow complexes are 
shown in Fig. |[ The source OH26.5+0.6 still shows an 
absorption spectrum in the 23 micron complex for both 
the amorphous and crystalline silicates, therefore we have 
neglected this source for the calculation of the mean out- 
flow spectrum. 

The complex is much less influenced by the presence of 
amorphous silicates and the crystalline silicate bands are 
more prominent both for the disk and outflow sources. 
The complex can be divided into 6 (blending) features 
at 22.4, 23.0, 23.7, 23.9, 24.5 and 25.0 /zm. These last 
two features form the 25 micron plateau. The 23.9 micron 
feature is only seen in the spectra of MWC 922, AC Her 
and OH26. 5-1-0.6 (see Table | and Paper I). In the latter 
object, it is present in absorption. 

The main difference between the mean complexes of 
the disk and the outflow sources is the relative strength 
of the 23.7 micron feature compared to the 23.0 micron 
feature, which is higher in the disk sources than in the 
outflow sources. We note that variations exist from source 
to source within each group. 

Another difference is that the spectra of disk sources 
show a strong 24.5 micron feature in the 25 micron 
plateaus, while the spectra of outffow sources do not. A 
third difference is the steep rise of the 23.0 micron fea- 
ture, present in most outffow sources, but not in the disk 
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Fig. 4. The normalised mean 28 micron complex spectrum 
of the disk and outflow sources. All sources are used, ex- 
cept Roberts 22 for the mean disk spectrum and AFGL 
4106, NML Cyg and IRC-hl0420 for the outflow sources, 
because of limited wavelength coverage. 

sources. Finally, we note that in almost all outflow sources 
there is a sharp drop around 22.5 /im, while this is not 
found or at least much less prominent in the disk sources. 

3.4- 28 micron complex (26.5 - 31.5 ^m) 

In Fig. ^, we show the mean 28 micron complexes for 
both the disk and outflow sources. For 4 sources, com- 
plete wavelength coverage was not available. These stars 
(Roberts 22, AFGL4106, NML Cyg and IRC-hl0420) have 
therefore not been taken into account in the mean disk and 
mean outflow spectra. This especially influences the mean 
outflow spectrum, since, apart from the 3 sources with- 
out a complete wavelength coverage, there is one source 
which goes from absorption into emission around these 
wavelengths and 2 sources are severely affected by noise. 
Therefore the mean outffow spectrum in this wavelength 
region should be treated with caution. This wavelength 
range is dominated by SWS band 3E (« 27.5 - 29.2 //m), 
whose ffux and relative spectral response calibration are 
somewhat more uncertain than for the other bands. Nev- 
ertheless, our large sample permits us to draw conclusions 
about the features present in this wavelength range. 

This complex is built up from features at 26.8, 27.6, 
28.2, 28.8, 29.6, 30.6 and 31.2 /im. The 29.6 /xm feature 
might be a blend itself, as indicated by the large variation 
in the peak position (Table |l|). However, one should be 
careful since there is an SWS band edge in the same region. 

The 29.6 and 30.6 micron features seem related to 
each other. A strong 30.6 micron features corresponds to 
a strong 29.6 micron feature and vice versa. 

3.5. 33 micron complex (31.5 - 37.5 ^m) 

In Fig. ^ the 33 micron complexes are shown. The shape 
of this complex bears similarities to the shape of the 
23 micron complex. Waters et al. (1996) noted the pres- 
ence of 3 main components: at 32.8 /im, at 33.7 /zm and 
a plateau extending up to 37 /im, hereafter called the 
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Fig. 5. The normalised mean 33 micron complex spectrum 
of the disk and outflow sources. 
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Fig. 6. The normaUsed mean 40 micron complex spectrum 
of the disk and outflow sources. 



35 micron plateau. However a closer look reveals much 
more structure. This is most clearly seen in the spectrum 
of NML Cyg (see Paper I). There is a weak feature at 
32.2 fim. The 33.7 micron peak is split into at least two 
features, peaking around 33.6 and 34.1 /xm. Since these 
two features are seen in both the AOTOl and the AOT06 
data of NML Cyg and in the spectra of HD179821 and 
AFGL4106 (see Paper I) we are convinced that these 
structures are real. The plateau shows some structure 
around 36 /im. It seems that it can be split into three 
features, at 34.9, 35.9 and 36.5 /im. These features are in 
most cases severely blended and just produce one broad 
feature (the 35 micron plateau). 

In the infrared bright outflow sources, the 32.97 mi- 
cron feature dominates, with IRC+10420 as the extreme 
case (see Paper I), while in the disk sources the 32.8 mi- 
cron feature is relatively strong, with AC Her as the ex- 
treme case (Paper I). However, the 32.97 micron feature 
coincides with a feature in the responsivity curve. It is 
most prominently seen in the high flux sources and in our 
sample the outflow sources are on average much brighter 
than our disk sources. For comparison we also checked 
the bright C-rich source CRL2688, which also shows the 
32.97 /im feature. It is also seen in the spectrum of other 
bright objects such as in Orion (Wright et al. 2000). This 
suggests a flux-dependent shape of the band 4 responsiv- 
ity curve, probably related to band 4 memory effects. We 
conclude that the 32.97 micron feature is instrumental. 

We note that the features seem sharper for the outflow 
sources than for the disk sources. The (post-)Red Super- 
giants, which have the highest luminosity and mass loss 
rates, seem to have the sharpest peaks. Why this is not 
the case for the other complexes is not clear yet. 

The detailed shape of the plateau, especially at the 
short wavelength side, is somewhat difficult to determine 
due to the presence of the strong bands between 32.5 and 
34.0 fim. The red edge of the plateau falls at 36.5 - 37 fim. 
This seems correlated with the sharpness of the strong 
bands in this complex. In the sources with the sharpest 
bands the plateau extends to 36.5 /im, while in most other 
sources the plateau extends to 37 /im. 



The apparent difference between the disk and outflow 
sources in the strength ratio of the 32.8 and 33.6 micron 
features is totally due to the (post-)Red Supergiants, since 
the low mass sources show about the same strength ratio 
as the disk sources. It might be a good discriminator be- 
tween (post-)AGB stars and (post-)Red Supergiants. 

3.6. 40 micron complex (38 - 45-5 fim) 

The mean 40 micron complexes (shown in Fig. ^ have 3 
prominent features, one at 40.5 /im and the often blended 
features at 43.0 and 43.8 //m. In addition to these main fea- 
tures, there is a 'plateau' visible from about 41 to 42.3 /im, 
hereafter called the 41 micron plateau, which is probably 
a blend of features, which could not be separated. Further- 
more, at the short wavelength side of the 40.5 /im feature, 
a shoulder is often seen around 39.7 /im and a weak fea- 
ture is present at 44.7 /im at the long wavelength side of 
the 43.8 micron feature . 

The main difference between the outflow- and the disk- 
sources is seen in the strength of the 43.0 and 43.8 micron 
features. In the outflow sources, the 43.0 micron feature is 
often stronger and at least equally strong to the 43.8 mi- 
cron feature, whereas in the disk-sources the 43.8 micron 
feature often dominates the 43.0 micron feature. 

Another difference between the disk and outflow 
sources is the blue edge of the 40.5 micron feature. In the 
outflow sources it has a more gentle slope and it meets the 
continuum at wavelengths shortwards of 39.0 /im, while in 
the disk sources the slope is much steeper and ends usually 
at wavelengths longwards of 39.0 /xm. 

3.7. 60 micron complex (50 - 72 fim) 

For this part of the spectrum, not all sources were thor- 
oughly reduced. Therefore, we will only mention the global 
trends and not go into details. The 60 micron complex 
(shown in Fig. ^ consists of a broad band peaking at 60 
- 65 /xm with a shoulder at 52 micron. In all disk sources, 
except for Roberts 22 (see Paper I), we find a prominent 
and sharp 69 micron feature, while this feature is much less 
prominent in the outflow sources. It is interesting to note 
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Fig. 7. The normalised mean 60 micron complex spectrum 
of the disk and outflow sources. The normalised mean 
spectrum for the disk sources is derived from NGC6537, 
NGC6302, MWC922, HD44f79 and Roberts 22, the nor- 
malised mean spectrum for the outflow sources is derived 
from all outflow source spectra. 

that Roberts 22 also shows similarities to outflow sources 
in other respects, i.e. low crystalline over amorphous ratio 
(see MYW). 

There are indications of a feature at 57 /im, but this 
has not been measured separately. The strong, narrow 
peak seen at 63. f 8 /im is due to the fine- structure transi- 
tion of [O I]. 

Disk and outflow sources differ in the peak position of 
the 60 - 65 micron band. In the outflow sources, it gen- 
erally peaks around 62 //m, while in the disk-sources it 
peaks at about 65 micron (with Roberts 22 as the excep- 
tion). The shift to longer wavelengths seems related to a 
decreasing 52 micron shoulder, and increasing strength of 
the 29.6 and 30.6 micron features. 




20 20.5 21 21.5 22 25.5 26 26.5 




wavelengtfi (/U,m) 

Fig. 8. An overview of the remaining features, the dashed 
line is the normalised mean spectrum for disk sources and 
the dotted line for outflow sources. There is no 90 micron 
feature detected in the outflow sources. 



3. 8. Remaining features 

In Fig. H, we show the mean spectra of the other features 
which were not attributed to one of the complexes. Some 
of these features tend to appear in most stars (C-rich and 
0-rich) and always in emission like the ones at 11.05, 13.6 
(with a shoulder at 13.8) and 14.2 micron and are there- 
fore believed to be (partly) instrumental. The strength of 
these features differs from source to source. However, it 
should also be noted that in some C-rich stars some of 
these features are seen up to 50% times the continuum 
level, which is far more than expected for an instrumental 
artifact. In these cases the features are attributed to PAHs 
(Tielens et al. 1999). For most of our sources however, we 
can exclude this possibility, since there is absolutely no 
evidence for PAHs at shorter wavelengths, which should 
be prominent if PAHs are present. 

Real features are present at 20.7, 21.5, 26.1, 47.7 (with 
a shoulder at 48.7) and 91 /im. 

The objects HD45677 and OH26. 5-1-0.6 both show a 
feature around 38 /im (see Paper I). The width of these 
two features is quite similar, but the peak positions differ 
about 0.4 /im. Because of the wavelength difference and 



because it is not seen in any of the other stars in our 
sample, makes us doubt its reality. 

4. Identification 

In this section, we will first discuss some of the laboratory 
measurements of crystalline silicates that may be present 
in circumstellar shells. Waters et al. (1996) identified most 
of these features with olivines (Mg2a;Fe(2_2a:)Si04) and py- 
roxenes (Mg^^Fei-jjSiOa), with 1 > x > 0. It turned out 
that the best match with t he o bserved features was found 
for X — 1 (see also Section 18); i.e. the minerals forsterite 
(Mg2Si04) and enstatite (MgSiOs). We will show the sim- 
ilarities and the differences between different laboratory 
data sets of these materials. The origin of these differences 
will be discussed. Finally, we will discuss the identification 
for the different features per complex. 

Laboratory measurements of crystalline silicates 

In order to understand and identify the different fea- 
tures, which were found in the infrared spectra of the dust 
around evolved stars, good laboratory spectra are neces- 
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Fig. 9. The mass absorption coefficients of forsterite derived from different laboratory measurements by, Koike et al. 
(1999; solid line), Jager et al. (1998; dotted line), Koike et al. (1993; dashed line). The Jager et al. data has been 
multiplied by a factor 2 to match the other measurements 



sary. Most laboratory spectra do agree on a global scale, 
but detailed comparison often reveals quite some differ- 
ences. In Fig. ^, [l^ and |ll|, we show different laboratory 
results for materials with the same bulk composition. It 
is clear that, although the average trends are very sim- 
ilar, there are significant differences in the details; e.g., 
strength ratios. 

The general trend in the 3 data sets of forsterite (taken 
from Jager et al. 1998, Koike et al. 1993 and Koike et 
al. 1999) in Fig. || is the same. Forsterite has two strong 
features in the 10 micron region, two prominent features 
in the 18 micron region, a strong band in the 23 micron 
complex and a strong band in the 33 micron complex. Fur- 
thermore, there is a weaker contribution in the 28 micron 
complex. A more detailed examination reveals quite some 
differences between the laboratory data sets. Peak posi- 
tions are not identical, this might be a grain shape effect, 
since the position (and shape) of the forsterite features 
are grain shape dependent. In particular, the difference in 
the peak position of the 23.7 micron feature for spherical 
particles and for a continuous distribution of ellipsoids is 
almost 1 fim. Therefore small differences between the lab- 
oratory spectra and the ISO spectra can be expected and 
do not lead automatically to a rejection of the proposed 
identifications. The Jager et al. data set of forsterite shows 
extra features at 9.3, 40.7 and 43.3 /zm. These are prob- 
ably due to enstatite inclusions. The measured absolute 
and relative strength of the features also differ. Note that 
we rescaled the Jager et al. data by a factor 2, in order 
to match the other studies. At the moment, the origin of 
this difference is unclear, but introduces substantial un- 
certainties in the derived relative abundances (paper III). 



For enstatite there are two crystallographic configu- 
rations, the rhombic structure (ortho-enstatite) and the 
monoclinic structure (clino-enstatite). On Earth, ortho- 
enstatite is most common, but enstatite found in mete- 
orites and produced via melting in laboratories is often 
of the clino-enstatite type. The structure of the enstatite 
grains is related to its thermal history, which makes them 
in principle a suitable candidate to explore the dust condi- 
tions when these grains were formed. Below 40 /im, there 
is not much difference between the spectra of clino- and 
ortho-enstatite. The main differences are found in the 40 
and 60 micron complexes. 

In Fig. |l0|, we plot the mass absorption coefficients 
of clino-enstatite derived from three different laboratory 
samples (Jager et al. 1998; Koike & Shibai 1998 and Koike 
et al. 1999). Again, the general trend is very similar: a lot 
of structure in the 10 micron complex, at least 7 differ- 
ent features can be found here, a strong feature peaking 
slightly longwards of 19 /zm and decreasing to longer wave- 
lengths with a lot of substructure. Enstatite has several 
modes in the 40 micron complex, where forsterite is rela- 
tively smooth. In this wavelength region, it also becomes 
possible to distinguish between clino- and ortho-enstatite, 
which is very difficult below 40 /im, and this spectral re- 
gion provides a prime probe for enstatite particles. Un- 
fortunately, the data sets do not agree in this wavelength 
region. The sample preparation of the Jager 1998 sam- 
ple is likely responsible for the absence of the 44.7 micron 
feature in their data (Jager, private communication). A 
careful examination of the ISO spectra with the labora- 
tory data convinced us that the Koike 1999 data set (Koike 
et al. 1999) best represents the features found in our ISO 
spectra. 
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Fig. 10. The mass absorption coefficient of clino-enstatite derived from different laboratory measurements by, Koike 
et al. (1999; solid line), Jager et al. (1998; dotted line), Koike & Shibai (1998; dashed line). This last data set only 
extends to 39 /xm. 



The spectrmn of ortho-enstatite (Fig. [Tl|; data from 
Jager et al. 1998, Koike & Shibai 1998 and Koike et al. 
1999) is a clear example why the purity of the laboratory 
samples always should be checked. The data from Jager 
et al. (1998) is contaminated by talc, producing various 
spurious features. Up to 40 /im, the ortho-enstatite is very 
similar to the clino- enstatite. The feature at 33.0 fim in 
the Koike 1993 data set (Koike et al. 1993) could not be 
reproduced by other measurements and is therefore not 
believed to be real. Therefore, the identification of ortho- 
enstatite to the 32.8 micron feature by Waters et al. (1996) 
is reconsidered. 

Diopside (MgCaSi206) is another silicate with possi- 
ble astrophysical relevance. This material shows peaks at 
19.5, 20.6, 25.1, 29.6, 32.1, 33.9, 40.1, 44.8 and 65.7 ^iul. In 
Fig. we show two laboratory measurements of this ma- 
terial (Data from Koike et al. 1999 and Jager private com- 
munication). Again differences are found between these 
two data sets, especially around 30 /im. 

These four examples show that one should be careful 
when applying a data set from the literature, especially if 
it has not been checked on chemical homogeneity, particle 
size and shape distribution. The homogeneity is not only 
important for the natural samples, which often contain 
inclusions of other materials (see e.g. Jager et al (1998)), 
but synthetic samples are also not always homogeneous. 

It is important to realize that differences in the forma- 
tion history of the minerals may lead to different struc- 
tures and IR spectral characteristics. Enstatite crystals 
which are formed from the vapour phase in a low density 
hot gas environment often produce structures in which one 
of the crystallographic axes is severely depressed (Bradley 
et al. 1983 and references therein). This is thought to be 



the situation in the outflows of mass losing stars. Therefore 
the absence or weakness of one of the laboratory bands in 
the stellar spectrum might be related to this phenomenon. 

In the following subsection, we will use the Koike 1999 
data set for the identifications, since that is thought to 
be the best and purest. However, one should realize that 
small wavelength shifts (2%) are found between the dif- 
ferent laboratory samples, all of them claimed to be pure 
forsterite. Small defects in the ordering of the atoms might 
cause these differences. It is likely that similar differences 
are also found for the grains in space. One should also re- 
alize that the laboratory data sets were measured at room 
temperature, while most dust grains are much colder. A 
lower dust temperature often shifts the wavelengths of the 
peaks to shorter wavelengths and strengthen the features 
(see e.g. Henning & Mutschke 1997, Bowey et al. 2000). 

4-2. 10 micron complex (1 - 13 iim) 

In Fig. |l^, we show the mean disk and outflow source spec- 
tra together with the Koike 1999 data sets. We subtracted 
a continuum from these data sets, in the same way as we 
did for the infrared spectra of our sources, in order to get a 
better comparison between the laboratory and ISO data. 
This continuum subtracted data set is also used in the 
other complexes. Note the similarities between the clino 
and ortho-enstatite samples. 

The broad feature in the 10 micron complex is thought 
to originate from amorphous silicates. This feature peaks 
at a relatively long wavelength (for both the average disk 
and outflow 10 micron complex) compared to the interstel- 
lar medium amorphous silicate absorption feature (which 
is at 9.7 /im). This can be an indication for the presence 
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Fig. 11. The mass absorption coefficient of ortlio-enstatite derived from different laboratory measurements by, Koike 
et al. (1999; solid line), Jager et al. (1998; dotted line), Koike & Shibai (1998; dashed fine). The Jager et al.(1998) 
data is contaminated with talc at 9.8, 14.9 and 21.5 and 25.4 /im. 
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Fig. 12. The mass absorption coefficient of diopside derived from different laboratory measurements by, Koike et al. 
(1999; solid line) and Jager (private communication). 



of large grains in the dust around evolved stars (see e.g. 
Molster et al. 1999b). In Fig. |l| we show this effect for 
spherical olivine grains (MgFeSi04; Dorschner et al. 1995). 
Note the shift in peak position and the broadening with 
increasing grain size. Besides grain size, the peak position 
and width of a feature can also be influenced by compo- 
sitional changes (see e.g. Tielens 1990; Dorschner et al. 
1995; Demyk et al. 1999). The redder peak position of the 
outflow sources is therefore probably not simply due to the 
presence of larger grains. NB, it is the disk sources which 



do show evidence for very large grains (see e.g. Molster et 
al. 1999a). 

In the mean disk spectrum, several sharp features are 
present which might be identified with crystalline silicates. 
The forsterite feature, found in the laboratory data at 
11.3 /im, is slightly red-shifted, while the other strong 
forsterite feature, at 10.0 /xm in the laboratory data, is 
blue-shifted. These shifts are within the differences seen 
between the different laboratory data sets. According to 
the laboratory data enstatite would peak between and 
around the forsterite peaks, causing a severe blend of the 
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Fig. 13. The normalised mean 10 micron complex spec- 
tra together with the forsterite (F; -I-2.3E4) and clino- 
enstatite (C; -t-5E3) and ortho-enstatite (O) mass ab- 
sorption coefficient spectra in this wavelength range. The 
forsterite and clino enstatite spectra are shifted. The gray 
lines indicate the features found in our spectra, with their 
spread in peak position observed in the ISO data indicated 
above the plot. The amorphous silicate peak has not been 
indicated. 



19 16 17 

wavelength (fJ-m) 



19 20 



Fig. 15. The normalised mean 18 micron complex spec- 
tra together with the forsterite (F; -I-1.3E4) and clino- 
enstatite (C; -f6E3) and ortho-enstatite (O) mass ab- 
sorption coefficient spectra in this wavelength range. The 
forsterite and clino-enstatite spectra are shifted. The gray 
lines indicate the features found in our spectra, with their 
spread in peak position observed in the ISO data indicated 
above the plot. The amorphous silicate peak has not been 
indicated. 
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Fig. 14. The Qabs for spherical olivine grains (MgFeSi04; 
Dorschner et al. 1995) for different grain sizes: 0.01 /im 
(solid line), O.lfim (dotted hne), 0.3/im (short dashed 
line), 1.0fj.m (long dashed Hne), S.Ofiui (dashed dotted 
line). The curves were scaled to get an equal strength 
around 10 fim and then offset from each other. 

features which hampers identification of the separate fea- 
tures. The strong band at 9.3 fim in the enstatite labora- 
tory data seems split into two components in the stellar 
spectrum. A hint of this splitting is already seen in the 
laboratory data. The 10.7 micron band in the ISO data 



coincides nicely with one of the strongest enstatite bands. 
In the mean disk spectrum, the 11.4 micron band is likely 
a blend of enstatite (at 11.7 fim) and the strong forsterite 
feature (at 11.3 /im). This blend might also explain the 
(apparent) red-shift of the forsterite band in our ISO-data. 

4-.S. 18 micron complex (15-20 ^m) 

The mean spectra of the 18 micron complex are shown in 
Fig. ^ together with the laboratory spectra of forsterite 
and both types of enstatite. Note that the difference be- 
tween clino and ortho-enstatite is only marginal, implying 
that these two forms cannot be distinguished in this wave- 
length region. 

Forsterite shows features around 16.3, 18.2 and 19.6 
and a shoulder at 19.0 /im (Jager et al. 1998; Koike et al. 
1999;), while enstatite shows features around 15.4, 17.6, 
18.5 and 19.5 /im (Jager et al. 1998; Koike et al. 1999). The 
laboratory data of these silicates are only slightly shifted 
from the peak positions found in our ISO-data. Because 
all expected peaks are found, not only in this wavelength 
range but also in the other complexes, the identification 
of enstatite and forsterite is unambiguous. However, the 
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Fig. 16. The normalised mean 23 micron complex spec- 
tra together with the forsterite (F; -I-1.1E4) and clino- 
enstatite (C; -t-7E3) and ortho-enstatite (O) mass ab- 
sorption coefficient spectra in this wavelength range. The 
forsterite and clino-enstatite spectra are shifted. The gray 
lines indicate the features found in our spectra, with their 
spread in peak position observed in the ISO data indicated 
above the plot. 

peak strength of the features is not as expected. Accord- 
ing to the laboratory data the peak at 19.5 ^m should 
have been by far the strongest feature in this complex. 
This is not observed and cannot be explained by temper- 
ature effects. Since there are also unidentified features at 
16.0 and 16.9 ^m it is likely that more dust components 
are present. These dust species might also contribute at 
other wavelengths in this complex and therefore change 
the apparent strength ratios of the different features. In 
particular, the 18.1 and 18.9 micron features likely have 
a component underlying the forsterite and enstatite fea- 
tures. 

The 16.5 micron feature is only seen in two stars 
(NGC6537 and HD44179, see Paper I), which both show 
PAH emission at the shorter wavelengths. There is a 
known PAH feature at about 16.4 /im (Van Kerckhoven 
et al. 2000), so we tentatively identify this feature with 
PAHs. 

4-4- ^3 micron complex (22 - 25.5 fim) 

The mean spectra of the 23 micron complex are shown in 
Fig. ^ together with the laboratory spectra of forsterite 
and both types of enstatite. Clino and ortho-enstatite 



show similar weak spectral structure, but the strength ra- 
tio of the peaks is reversed. 

Both forsterite and enstatite have solid state bands 
in this wavelength range. These two components can ac- 
count for some of the observed structure. The 23.7 micron 
feature has been attributed to forsterite, while the 23.0 
and 24.5 micron features coincide with enstatite features. 
However, since the strength of these two bands are not cor- 
related (in the outflow sources the 24.5 micron feature is 
weak and the 23.0 micron relatively strong) it is expected 
that, at least one, of these two features has a contribution 
from another dust species. 

The 22.4 and 25.0 micron features still lack identifica- 
tion, they seem not to be correlated with the other features 
in this complex. In most outflow sources we see a drop 
around 22.5 fiui. Since absorption at these wavelengths 
is not expected (for most sources), it is either caused by 
emission features at both sides, or due to responsivity. 
Since it is not seen in the disk sources and the two emis- 
sion features are found in several sources, we attribute the 
22.5 iim drop to the narrowness of the crystalline silicate 
features in the outflow sources. 

4-5. 28 micron complex (26.5 ~ 31.5 fj,m) 

The mean spectra of the 28 micron complex are shown in 
Fig. ^ together with the laboratory spectra of forsterite 
and both types of enstatite. In this complex the spectra 
of ortho and clino-enstatite start to deviate slightly from 
each other. 

Laboratory measurements show that forsterite has 
two features in this wavelength range: a strong band at 
27.5 fim and a weaker (w 1/3 of the 27.5 /xm) feature 
at 31.2 /im. This second feature seems to disappear if Fe 
is included in the lattice structure (see Jager et al. 1998; 
Koike et al. 1993). The spectrum of NGC6537 (see Pa- 
per I) is remarkable in this respect, since it does show 
the typical forsterite features and even a 69.0 micron fea- 
ture (indicating that the forsterite is very magnesium rich, 
(Sec ^.7| and Jager et al. 1998), while there is no indication 
of a 31.2 micron feature. It is even more remarkable con- 
sidering the fact that the spectrum of NGC6537 is very 
similar to that of NGC6302 (Paper I) which shows the 
most prominent 31.2 micron feature in our sample. This 
raises some doubt about the forsterite identification of the 
31.2 micron feature in the ISO spectra and we leave this 
to future research. 

The 28.2 micron feature, often blended with the 
27.5 micron feature, is most likely due to enstatite, al- 
though it is not clear why there is no evidence for an 
almost equally strong 26.5 micron feature, which is seen 
in the laboratory transmission spectra (Jager et al. 1998; 
Koike et al. 1999). The weak 26.9 and 28.7 micron features 
as well as the relatively strong 29.6 and 30.6 micron fea- 
tures are not well described by forsterite and/or enstatite. 
For the latter two features this is very clear in the fit to 
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Fig. 17. The normalised mean 28 micron complex spec- 
tra together with the forsterite (F; -I-8.5E3) and chno- 
enstatite (C; -I-3.5E3) and ortho-enstatite (O) mass ab- 
sorption coefficient spectra in this wavelength range. The 
forsterite and clino-enstatite spectra are shifted. The gray 
lines indicate the features found in our spectra, with their 
spread in peak position observed in the ISO data indicated 
above the plot. 
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Fig. 18. The normalised mean 33 micron complex spectra 
together with the forsterite (F; -I-6E3) and clino-enstatite 
(C; -I-3E3) and ortho-enstatite (O) mass absorption co- 
efficient spectra in this wavelength range. The forsterite 
and clino-enstatite spectra are shifted. The gray lines in- 
dicate the features found in our spectra, with their spread 
in peak position observed in the ISO data indicated above 
the plot. 



the spectra of NGC6537 and NGC6302 (Paper III). The 
strength of 29.6 and 30.6 micron features are correlated. 
This would imply that they are originating from the same 
dust species. No identification has been made yet. 

4-6. 33 micron complex (31.5 - 37.5 fim) 

The mean spectra of the 33 micron complex are shown in 
Fig. ^ together with the laboratory spectra of forsterite 
and both types of enstatite. There is a small shift between 
the clino and ortho-enstatite features, which is most clear 
for their peak around 36 fim. 

Both forsterite and enstatite have features in this 
wavelength range. The band at 33.6 fim is identified with 
the strong resonance peak of forsterite. The strongest en- 
statite peaks are found around 36 ^m (35.3 fim for clino- 
enstatite and 36.1 ^m for ortho-enstatite). They might 
be the relatively warm carriers of the cold circumstel- 
lar 34.9 and 35.9 micron features. They are likely the 
main contributors to the plateau (see also the fits in pa- 
per III). The 36.5 micron peak also falls into the plateau 
regime. Forsterite does show a peak at these wavelengths, 
which would suggest that the long wavelength end of the 
plateau should correlate with the strength of the 33.6 mi- 



cron forsterite peak. However, this correlation is not found 
in our sample. Also, the strength of the laboratory 36.5 mi- 
cron feature is too weak to explain the strength of the 
36.5 micron feature seen in some sources, e.g. MWC922 
(Paper I). The fits to the plateau in Paper III make clear 
that we are also still missing a component that peaks 
around 34.5 fim. However we were not able to extract 
that feature from our spectra, which indicates that is is 
probably very blended with other features. 

If the red-edge of the plateau is due to ortho-enstatite, 
then the difference in the extension of the plateau (36.5 /xm 
for the outfiow objects and 37 /.tm for the disk sources) 
could be explained by a larger percentage of clino versus 
ortho enstatite. However, this is in contrast with what has 
been found in the analysis of the 40 micron complex (see 
Sec ^). So, we cannot be conclusive about this. 

Laboratory spectra of enstatite do show peaks at 33.0 
and 34.0 /xm, but these are much weaker than the en- 
statite features seen in the plateau. This is probably not 
a big problem for the 34.1 micron feature. Although the 
measured 34.1 micron feature is narrower than seen in the 
laboratory, this is the case for more features and might be 
a temperature effect. The 34.1 micron feature is a typical 
example of a feature that often does not show a sharp peak 
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Fig. 19. The normalised mean 40 micron complex spec- 
tra together with the forsterite (F) and clino-enstatite (C; 
-I-3.5E3) and ortho-enstatite (O; -I-1.2E3) mass absorption 
coefficient spectra in this wavelength range. The ortho- 
and clino-enstatite spectra are shifted. The gray lines in- 
dicate the features found in our spectra, with their spread 
in peak position observed in the ISO data indicated above 
the plot. 



but seems necessary to add to fit the spectrum. Since there 
is no correlation found with the 33.6 or 34.9 micron fea- 
ture it is not expected to be related to the non-Gaussian 
shape of these features. We note that the 32.97 micron 
feature in high flux sources is due to incorrect responsiv- 
ity curves (see Section 3.5). . The 32.2, and 32.8 micron 
features still lack solid identification. 



(see Paper I). This may indicate that two different dust 
components are responsible for the emission in the 40.5 /im 
range. It should be noted that these two components are 
not seen separately, which indicates that, if present, they 
severely blend. Although it cannot be excluded at the mo- 
ment that it is just one solid state, which can shift quite 
dramatically, we have other indications for a second com- 
ponent in the 40.5 micron feature. The fit to the spectra 
of HD179821 and AFGL4106 in Paper III demonstrate 
that in these cases another dust component is necessary 
to explain the strength of the features. 

The 43.8 micron feature can probably be explained 
with ortho-enstatite. The strength ratio between the 40.5 
and 44.7 micron features might therefore give an indica- 
tion of the clino-enstatite over ortho-enstatite ratio. How- 
ever, one would expect that this difference is also seen at 
the long wavelength edge of the 35 micron plateau and 
there is no correlation found in that spectral regime. This 
indicates that other dust species influence the features at- 
tributed to enstatite features, as was already indicated for 
the 40.5 micron feature. 

The 43.0 micron feature is attributed to crystalline 
H2 0-ice and clino-enstatite. It should be noted that the 
wavelength of the crystalline H20-ice feature can shift sig- 
nificantly with temperature, from 42.7 at 10 K to 44.7 at 
140 K (Smith et al, 1994). Our peak positions indicate 
that the average ice temperature is around 40 K. This is 
a rather low temperature, although similar temperatures 
have been found for the Il20-ice in AFGL4106 (Molster et 
al 1999b) and HD161796 (Hoogzaad et al. in prep.), based 
on a comparison of the strength and position of the 43 and 
60 micron bands of H2 0-ice. The relative strength of the 
43.0 micron feature in the outflow sources compared to 
the disk sources suggests that crystalline Il2 0-ice is more 
abundant in the outflow sources than in the disk sources. 

The 39.7 micron feature lacks identification, but there 
are indications that its strength shows an anti-correlation 
with the strength of the 44.7 micron feature. We therefore 
speculate that this feature is due to a reaction product of 
clino-enstatite and another (dust) species. 



4-.7. 40 micron complex (38 - 45-5 fim) 

In Fig. ^9] we show the mean spectra and the laboratory 
spectra of forsterite and both types of enstatite. Forsterite 
shows no features in this wavelength region. The clino- 
enstatite feature at 44.9 /im, which we identify with the 
44.7 micron feature, is the first feature which is completely 
different from the ortho-enstatite features and therefore a 
good benchmark for the clino- to ortho-enstatite ratio. 

The 40.5 micron feature is probably caused by en- 
statite, since laboratory spectra of both ortho- and clino- 
enstatite show a peak at this position. The 40.5 micron 
feature shows large source to source variations in peak po- 
sition for such a well defined feature, the extremes being 
IRAS09425-6040 (40.35 urn) and HD179821 (40.80 ^ira) 



4-8. 60 micron complex (50 - 72 fj,m) 

Fig. ^ shows the mean spectra of the 60 micron complex 
and the corresponding spectra of forsterite and enstatite. 
We note that these silicates only show sharp features in 
their lab spectra, while this complex is dominated by a 
broad feature. 

The 62 micron feature and the shoulder at 52 /im, 
which is hardly seen in the disk sources and much more 
prominent in the outflow sources, is identified with crys- 
talline Il2 0-ice. In the disk sources the 62 micron band 
is red-shifted compared to the outflow sources. Labora- 
tory data of crystalline Il2 0-ice show only a small shift 
in peak position of the 60 micron band with temperature. 
This shift is too small to explain the large differences in 
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Fig. 20. The normalised mean 60 micron complex spec- 
tra together with the forsterite (F; -I-2.7E3) and chno- 
enstatite (C; +1.6E3) and ortho-enstatite (O) mass ab- 
sorption coefficient spectra in this wavelength range. The 
forsterite and chno-enstatite spectra are shifted. The gray 
lines indicate the features found in our spectra, with their 
spread in peak position observed in the ISO data indicated 
above the plot. 



peak position between the disk and the outflow sources. 
This indicates that the 62 micron feature actually contains 
at least two solid state features: a 60 micron feature at- 
tributed to crystalline water ice and a 65 micron feature. 
The shift to shorter wavelengths in the outflow sources im- 
plies that the relative abundance of crystalline water ice 
is higher in the outflow sources, a conclusion which was 
also drawn on the basis of the 40 micron complex. The 
65 micron feature is well fitted with diopside (Koike et al. 
2000). Enstatite also shows some peaks in the 60-70 /im 
wavelength range, but since they are quite sharp in the 
laboratory data (Koike et al. 1999) it is not likely that 
enstatite is the main contributor to the 65 micron band. 

The 69.0 micron feature has been identified with 
forsterite. This peak shifts significantly by adding fayalite 
to the solid solution series (already to 73 ^m for a com- 
position of 4% fayalite and 96% forsterite, see Jager et 
al. 1998). This feature therefore indicates that the olivine 
dust particles are very pure forsterite with hardly any Fe. 
However, it should be noted that the laboratory peak po- 
sition of pure forsterite at room temperature is 69.7 /xm, 
while we observe this feature around 69.0 /im. Recent lab- 
oratory measurements at low temperatures (Bowey et al. 
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Fig. 21. A comparison of the features around 14 /im 
in the continuum subtracted spectra of OH26.5+0.6 and 
AFGL4106, together with the responsivity proffie and the 
mirrored (in the wavelength-axis) mass absorption coeffi- 
cient spectrum of enstatite. Three dips in the responsivity 
are indicated by dashed lines. 

2000,2001; Chihara et al. 2001) indicate a shift to shorter 
wavelengths with decreasing temperature; e.g. at 4 K this 
feature narrows significantly and peaks at 68.8 /im. This 
would provide a natural solution to the discrepancy be- 
tween observed and laboratory data taken at room tem- 
perature, see also Paper IIL 

4-9. Remaining features 

It was noted that the features around 14 fim might be 
(partially) instrumental. The best indication for this state- 
ment would seem OH26.5-1-0.6 which shows these spurious 
features also in emission while all other features in this 
wavelength range are in absorption. In Fig. ^ we show 
its spectrum in this wavelength region, together with the 
spectrum of AFGL4106. 

We would like to point out that the spectrum of 
OII26.5+0.6 could be very well fitted by the 'absorption' 
spectrum of enstatite. However the features around 14 /im 
were also found in stars which only show emission like 
AFGL4106. In these cases the enstatite spectrum alone 
can not explain the structure anymore, and it looks more 
like a responsivity problem, also because it shows quite 
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some similarities with the responsivity spectrum. Caution 
should therefore be taken for these features, since the are 
likely influenced by both enstatite and responsivity prob- 
lems. 

The 20.7 micron feature is seen in the forsterite and 
enstatite data, but not as strong as it is seen in some 
sources. Silica (Si02) does show a very strong resonance at 
this wavelength and might therefore be a good candidate. 
If silica is present it can be an explanation for the 9.1, 15.9 
and 26.1 micron features (Fabian et al. 2000). However, 
one would also expect to see a feature at 12.6 micron which 
is not the case. Therefore, a solid identification for this 
feature is still lacking. 

The 21.5 micron feature is also seen as a weak feature 
in the enstatite and forsterite laboratory spectra. In some 
sources it is stronger than expected from these laboratory 
spectra. Therefore, it is likely that another material con- 
tributes to this feature. While the 20.7 micron feature is 
always present when the 21.5 micron feature is detected, 
this is not true the other way around. Therefore, they must 
originate from different dust species. Also, the strength ra- 
tio is not always the same. 

Forsterite has a feature at 26.1 /im (Jager et al. 
1998). We, therefore, identify the 26.1 micron feature with 
forsterite, although a contribution from silica cannot be 
excluded. 

The blend of the 47.7 and 48.6 micron features has 
not yet been identified; enstatite nor forsterite have fea- 
tures at these wavelengths. Ferrarotti et al. (2000) iden- 
tified the blend of the 47.7 and 48.6 micron band with 
FeSi. The 44.7 micron features might also be explained by 
this material. However, their fit is not completely satis- 
factory and the FeSi features around 30 fxm are also not 
well matched with features in our ISO spectra. We wonder 
if shape effects shifts (and reshapes) the bands as is seen 
for forsterite. This has not been tested. This identification 
was based on a condensation sequence in stellar outflows 
where free carbon and oxygen atomes are both underabun- 
dant. This can be either because the C/0 ratio is close 
to 1 and all carbon and oxygen is locked up into CO, or 
in stars which expose their layers, which formerly have 
burnt hydrogen by the CNO-cycle. However, we observe 
the 48 micron feature in several stars, some of which are 
definitely 0-rich (e.g. OH26.5+0.6). The proposed identi- 
fication of the 48 micron blend with FeSi is an interesting 
identification, but not convincingly confirmed yet. Other 
(hydrous) silicates often have peaks around 48 /im. How- 
ever, no suitable candidate has yet been identified. There- 
fore, the identification of the 47.7 and 48.6 micron features 
remain open. 

The longest wavelength feature which has been found 
in various stars is the 91 micron feature. Malfait et al 
(1998) suggested hydrous silicates; these materials often 
show a broad feature in this wavelength range, but a sat- 
isfying fit has not yet been achieved. It is also difficult 



to explain how hydrous silicates would form in these out- 
flows, since liquid water seems necessary. 

In 89 Her and HD44179 a broad (full width zero inten- 
sity « 15 /im) solid state band which peaks around 33 /im, 
hereafter called the 33 micron band, is found. It has not 
been included in Table |[ because its shape and appear- 
ance are very uncertain. It shows a relatively sharp rise on 
the blue side at about 27 /xm while it ends with a gentle 
slope on the red side at about 42 /zm. In the amorphous 
silicate data presented by Ossenkopf et al. (1992) a broad 
feature seems present around these wavelengths, but the 
limited wavelength resolution of this data set at these 
wavelengths prevents a firm identification. The strength 
of this feature in our ISO data is likely too strong to be 
solely explained by the feature in the Ossenkopf data set. 
Therefore, other materials are likely to contribute to this 
feature. The peak position and width bear some similari- 
ties with the 30 micron feature in C-rich sources, where it 
is attributed to MgS. Similar features are also seen in the 
direction of the galactic center and in the dust of rj Car. 
Note that the shape of the 28 micron complex seems in- 
fluenced by this feature. 

5. Discussion and conclusions 

In this paper we have combined the spectra of 17 stars 
with oxygen-rich circumstellar dust environments. These 
combined spectra show a wealth of dust features, which 
now are studied in a systematic way. Most of the features 
cluster into one of seven complexes (at 10, 18, 23, 28, 33, 
40 or 60 /im) . Each of these complexes is a convolution of 
a limited set of components, which vary independently. 

Our conclusion from this work is that there is a clear 
distinction between the disk and outflow sources. Not 
only in the strength of the features and therefore in 
abundance of the crystalline silicates, which was the orig- 
inal natural separator profile between these two groups 
of objects, but also in the characteristics of the complexes. 

10 + 18 micron complex: 

The difference in shape, a smooth profile for the outflow 
sources and narrow features for the disk sources, is mainly 
due to an abundance difference between the crystalline 
and amorphous silicates in the disk and outflow sources 
and is also directly connected to the strength of the 
features. The outflow sources have a lower abundance of 
crystalline silicates. 

23 micron complex: 

The strength of the 23.0 micron feature with respect 
to the strength of the 23.7 micron feature is higher in 
the outflow sources. The 24.5 micron feature is more 
prominent in the disk sources. In Paper HI we will 
demonstrate that in the (massive) outflow sources, the 
ratio of enstatite over forsterite is higher. This would 
imply that the 23.0 micron feature results indeed from 
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enstatite and that although enstatite has a (weak) feature 
around 24.5 fiia, another material is (mainly) responsible 
for the 24.5 micron feature. 

28 micron complex: 

The 28 micron complex differs significantly from source to 
source and several outflow sources are affected by a lack 
of data. Therefore, no clear difference could be pointed 
out. 

33 micron complex: 

The features appear sharper in the outflow sources than 
in the disk sources. This might have to do with the 
formation history of the crystalline silicates (see Paper 
III). The (post-)Red Supergiants separate out by almost 
equally strong 32.8 and 33.6 micron features, even after 
subtracting the spurious 32.97 micron feature. This is not 
a characteristic of all outflow sources. 

4-0 micron complex: 

For the 40 micron complex the main difference is found 
in the strength ratio of the 43.0 and 43.8 micron features. 
In the disk sources, the 43.8 micron feature is more 
prominent and in the outflow sources the 43.0 micron 
feature. This has to do with the high abundance of 
crystalline Il20-ice in the outflow sources. 

60 micron complex: 

The main difference in this complex is found in the peak 
position of the 62 micron feature. This feature is probably 
a blend of crystalline Il2 0-ice, peaking around 60 fim 
(the outflow sources), and diopside plus possibly other 
crystalline silicates, peaking around 65 /im (the disk 
sources). Also, the 69.0 micron feature is much stronger 
in the disk sources than in the outflow sources, which is 
probably an abundance effect. 

It is likely that the different history of these dust 
grains are responsible for the spectral differences. MYW 
suggested that low temperature crystallization has taken 
place in the disk sources. In contrast, the crystallization 
of silicates in outflow sources is likely to have taken place 
close to the star; i.e., high temperature crystallization. 

Most features can be explained with forsterite and en- 
statite. It is clear that several other features lack a proper 
identification, and also the strength ratios are not always 
correct, suggesting the presence of more contributing ma- 
terials. Diopside has been mentioned but the strength of 
several features attributed to diopside suggests that it is 
only a minor component. 

The growth of enstatite crystals can occur in a pref- 
erential direction of one of the crystallographic axes. In 
chondritic porous interplanetary dust Particles (IDPs), 
which are believed to be the most pristine material of 
our solar system, enstatite has been found in which one 
or 2 directions of the crystallographic axes were severely 



depressed (Bradley et al, 1983). These enstatite crystals 
are likely formed from the vapour phase in a low density 
hot gas environment (Bradley et al, 1983 and references 
therein) , in conditions which are similar to the conditions 
in the outflows of stars. It is therefore possible that one or 
two of the crystallographic axes ([010] and [001] according 
to what has been seen in IDPs) are depressed. This will 
change the emission spectrum of this material in such a 
way that certain feature will be weak, corresponding to the 
depressed crystallographic axis. This scenario might help 
to understand why the relative strength of some features 
differ from the lab spectra. 

Another abundant species is II2O ice. Kouchi and 
Kuroda (1990) demonstrate that below 70 K the ice 
rapidly turns amorphous under the bombardment of UV- 
photons. Calculations based on their results and reason- 
able assumptions for the UV flux of the central stars would 
imply that all the crystalline water ice turns amorphous 
within a couple of days. Apparently the ice is efficiently 
shielded from the UV radiation of these stars. Also the low 
temperature of the ice is an indication for the shielding of 
the crystalline water-ice. If the material is very clumpy as 
is seen in e.g. the Helix Nebula (first seen by Baade and re- 
ported by Verontsov-Velyaminov (1968), it might provide 
a suitable environment to survive the UV radiation. 

In the outflow sources the abundance of crystalline wa- 
ter ice seems much higher than in the disk sources. Assum- 
ing that initially the same amount of H2O was formed, this 
difference points to a scenario in which crystalline water 
ice is destroyed on time scales longer than the outflow time 
scales. In the outflow of stars, density enhancements has 
been found, where molecules and dust are protected from 
the harsh UV radiation. This might explain the crystalline 
water-ice content even for the cold dust. Also inside disks 
particles are protected against the UV radiation, but the 
turbulence in disks likely causes that most dust particles 
will ultimately be exposed to the destructive UV radia- 
tion, long before the disk has been blown away. 

With about 80% of the features identified, see Ta- 
bic |l| for the identifications, the age of astromineralogy 
has now really started. Still, about 20% of the features 
lack a proper identification. There is also ample evidence 
that for a significant amount of the 'identified' features 
more dust components are necessary to explain the rela- 
tive strength of the infrared bands. New laboratory mea- 
surements of cosmic relevant dust species up to at least 
100 ^m are required to identify these features. We can 
now start to fully exploit the historical information of the 
conditions the dust experienced, by studying the composi- 
tion and chemical structure of the grains. The differences 
found between the different laboratory datasets, are a sig- 
nificant obstacle in the quantitative interpretation of the 
spectra. The properties of the sample (e.g. the stoichiome- 
try, size of the particles, number and size of the individual 
crystals in a particle) should all be well checked before and 
after the measurents of the optical constants. We would 
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also rccomcnd that the different astrophysical laborato- 
ries exchange their samples in order to have a better un- 
derstanding of the possible sources of differences in their 
datasets. 
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